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The uncertainties of the mean diameters of 
the nominal 1.0 |im SRM® 1690 poly- 
styrene spheres and of the nominal 0.3 |im 
SRM 1691 polystyrene spheres are 
recomputed using the current NIST 
Guidelines for computing uncertainty. The 
revised expanded uncertainty (approxi- 
mately 95 % confidence level) for SRM 
1690 polystyrene spheres is equal to 0.005 
|im compared to previous value of 0.008 
|im. The revised expanded uncertainty for 
SRM® 1691 is equal to 0.004 |xm com- 
pared to the previous value of 0.007 |xm. 
The major cause of the reduction in the 
uncertainty for the 1 .0 |xm spheres is from 
a decrease in the recomputed uncertainty 
of the refractive index of the polystyrene 



spheres. The 1 .0 |xm spheres were used in 
calibrating the electron microscope used to 
size the 0.3 |xm spheres, and the reduction 
in the uncertainty of 1.0 |xm SRM uncer- 
tainty was the biggest factor in the 
decrease in the uncertainty of the 0.3 |xm 
spheres. 
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1. Introduction 

The NIST Standard Reference Material 1690 (NIST 
SRM® 1690) consists of a nearly monosize suspension 
of 0.895 |J,m polystyrene spheres in water at a mass 
fraction of approximately 0.5 %. This standard is used 
in the certification of secondary standards and also used 
directly in the calibration of electron microscopes, of 
scanning surface inspection systems in the semicon- 
ductor industry, and of other particle sizing instruments 
when the most accurate sizing standards are needed. 

The certification of SRM® 1690 was based on the 
measurement of light scattering intensity versus scat- 
tering angle for a diluted suspension of the polystyrene 
spheres. Key features of the experiment [1] were the 
use of an intensity stabilized laser, an accurately 
indexed rotary table, and photon counting detection. 
The particle diameter was determined from a nonlinear 



least squares fit of the predicted scattering based on 
Mie theory and the measured data. 

The uncertainty analyses used in this study [1] is not 
consistent with the current NIST policy [2,3] for report- 
ing measurement uncertainty and result in overesti- 
mates of those uncertainties. The focus of this note is 
the recalculation of the uncertainty per the 1994 
Guidelines [2]. There is also an updated analysis of the 
refractive index component of the uncertainty. 

NIST SRM® 1690 was also used in the measurement 
of SRM® 1691 (0.269 |J,m diameter polystyrene 
spheres). The effect of the change in the uncertainty in 
SRM® 1690 on the uncertainty for SRM® 1691 is deter- 
mined. The uncertainty in the SRM® 1691 is also 
recomputed per the 1994 guidelines, since the method- 
ology used previously [4] is not consistent with current 
NIST policy 
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2. Calculation of Expanded Uncertainty 

The old procedure for computing the uncertainty is 
briefly reviewed and then the revised uncertainty analy- 
sis is presented based upon the 1994 Guidelines. 

2.1 Old Method 

The total uncertainty, {/^(oid), was computed by 
adding the random error, R, and the sum of the absolute 
values of the systematic errors. Kg/- 



U. 



T(old) 



^^ + Il 



(1) 



The random component of the uncertainty, R, was com- 
puted as the product of a coverage factor, k, for a 95 % 
confidence level times the uncertainty of the mean for 
10 repeat measurements of the mean, u,. 
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The quantity D„,- is the average diameter of the i-th 
sample, and D is the average of the 10 samples. The 
computed value of m^ equals 0.000229 |xm. The cover- 
age factor k for 9 degrees of freedom based on 
Student's t-distribution for "about 95 %" confidence 
interval is 2.32. Thus the value of i? is given by the fol- 
lowing: 



R = ku^ = (2.32)(0.000229) =0.00053 i m. 



(3) 



The systematic uncertainties are related to the parti- 
cle properties and the optical system. The uncertainties 



are expressed in terms of the effect on the particle 
diameter and the values are given in Table 1. The parti- 
cle related uncertainties include the refractive index of 
the spheres, Mbd the presence of about 1 % agglomerat- 
ed doublets, Ug2, and multiple scattering from the parti- 
cle suspension, Mgs- The optical system related uncer- 
tainties include the reflection from the glass cell, Ug^, 
the finite acceptance angle of the detector of about ±1°, 
Mb5, and the slight optical misalignment at zero angle, 
Mb6- These uncertainties, which were referred to as sys- 
tematic uncertainties in Mulholland et al. [1], are now 
classified as Type B uncertainties. These estimates are 
based on scientific judgment rather than based on sta- 
tistical methods, as is done for Type A uncertainties. 

The total uncertainty, {/^(oid) = 0.0074, is computed 
from Eqs. (1), (3), and the sum of the systematic uncer- 
tainties (see Table 1). 

2.2 New Method 

In 1994 the method for reporting uncertainties at 
NIST was unified [2] and aligned with the ISO Guide 
to the Expression of Uncertainty [3]. In this approach 
each component of uncertainty of a measurement result 
is represented by an estimated standard deviation, 
termed standard uncertainty with symbol u. There are 
two types of standard uncertainty. The first is comput- 
ed by statistical means such as the standard deviation of 
the mean of several repeat measurements and is termed 
a Type A standard uncertainty. The second is often 
based on scientific judgment using all the relevant 
information available and is termed Type B standard 
uncertainty. 



Table 1. Type B (systematic) Uncertainties for Measurement of SRM 1690 



Type B uncertainties 


Symbol 


Original value, 


|jjn 


Corrected value, |im 


Refractive index 


"Bl 


0.0030 




0.0020 


Particle doublets 


"B2 


0.0010 




0.0010" 


Multiple scattering 


"B3 


0.0010 




0.0010" 


Cell reflection 


"B4 


0.0010 




0.0006 


Finite acceptance angle 


"B5 


0.0005 




0.0003 


Optical misalignment 


"B6 


0.0004 




0.0004" 


Il"B.I 




0.0069 






Ik.)' 






0.0035 




0.0026 


" These val 


ues are now 


considered to correspond to standard uncertainties. 





The symbol 8i, used in Mulholland et al. [1], has been replaced with 

"Bi- 
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In the case of NIST SRM® 1690, the Type A standard 
uncertainty is the standard deviation of the mean of 10 
repeat measurements of the number mean diameter, 
which is defined as u, in Eq. (2) and found to have a 
value of 0.000229 nm. 

The Type B standard uncertainties for NIST SRM® 
1 690 consist of six components of uncertainty given in 
Table 1 [1]. 

Following the NIST Guidelines, the combined 
uncertainty, u„ is computed as the root-sum-of-squares 
of the Type A uncertainties and the Type B uncertain- 
ties. The basis of this approximation is that provided 
the variables are independent, the variance of a sum of 
independent variables is equal to the sum of the vari- 
ances. 



u^ + 



X"i 



= 0.0035. 



(4) 



The expanded uncertainty, U, which defines an inter- 
val having a level of confidence of about 95 % 
(95.4 %), is computed as U= ku^. The quantity k is the 
coverage factor and its value is dependent on the num- 
ber of degrees of freedom for u^. In the limit of infinite 
degrees of freedom, the value of k is 2.0. For a finite 
number of degrees of freedom, k is estimated as the t- 
factor from the Student's t-distribution based on the 
number of degrees of freedom and about a 95 % confi- 
dence interval. For a combined uncertainty arising from 
several components each with degrees of freedom v„ 
the effective number of degrees of freedom, v^s, is esti- 
mated using the Welch-Satterthwaite formula [2]: 



(5) 



4 4 

c, u. 



It is assumed that the number of degrees of freedom for 
each of the Type B terms in Eq. (4) is infinity so that the 
only term in the sum is the Type A uncertainty given by 
Eq. (2). In this case the sensitivity factor c, is unity, the 
term m, = u, = 0.000229 |J,m, and the degrees of free- 
dom, v„ is 9 (the number of repeat measurements minus 
one). The resulting value of V^g from Eq. (5) is 
4.8 X 10^ The corresponding value for k = 2.00. Given 
k, the value of the expanded uncertainty, U, is comput- 
ed as 0.0070 iim. 



2.3 New Method With Revised Type B 
Uncertainty Estimates 

Several of the systematic uncertainties estimated by 
Mulholland et al. [1] are over estimates. The term "at 
mosf ' is used in describing both the uncertainty associ- 
ated with the reflections from the glass cell, u^^, and the 
finite acceptance angle of the detector, u^^. To convert 
these estimates to Type B standard uncertainties, we 
treat each of these quantities as having equal probabil- 
ity over the respective ranges of ± 0.001 |J,m for Mb4 and 
± 0.005 |J,m for Mgj. For this rectangular probability dis- 
tribution, the standard deviation is u^/{3y'^. So both of 
these uncertainties are reduced to 0.58 times their pre- 
vious values, which corresponds to 0.0006 |J,m for u^^ 
and 0.0003 |J,m for Mgs (See corrected values in Table 

1-) 

The previous estimate of refractive index uncertain- 
ty corresponded to the range in the reported values 
from five studies [1]. This provides an over-estimate 
and the estimate could be revised by the method used 
above for u^^ and Mgs- Instead, we compute the uncer- 
tainty based on the single particle refractive index 
measurements by Marx and Mulholland [5] for the 
SRM® 1690 particles. This is a more accurate approach 
because two of the other four studies involved particles 
at least a factor of three smaller than the SRM® and the 
other two studies, which used a method similar to [5], 
did not include a quantitative uncertainty analysis. The 
measurement of the refractive index was based on 
measuring the light scattering versus angle from 30° to 
1 60° for a single, levitated SRM® sphere. The refractive 
index and particle size were determined from best fits 
of Mie theory predictions to the scattering data for the 
incident laser polarization direction both parallel to the 
scattering plane and perpendicular to the scattering 
plane [5]. The best fit was based on the maximum in the 
harmonic mean, l/Q^ + I/Q2, of the results for the two 
polarization directions. The quantities Q^ and Q2 are the 
sums-of-squares of the differences between the meas- 
ured and predicted scattering for the laser polarization 
direction parallel and perpendicular to the scattering 
plane. The resulting mean and standard deviation of the 
mean for measurements on eight separate particles is 
1.6121 ±0.0013. 

There are two sources of Type B uncertainty for the 
refractive index: uncertainty in the angle, Ug, and in the 
polarization direction, Up. There was a slight drift in the 
encoder angle readout of 0.08° over the time that the 
measurements were made. Including this drift in the 
numerical simulation of the light scattering, it was 
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found that the change in the refractive index was 0.005. 
We use this value as our estimate of Ug. It was found 
that for each polarization direction selected, there was 
about 0.5 % of that intensity of light with the orthogo- 
nal polarization direction. From numerical simulations, 
this effect was found to change the refractive index by 
0.0032 (0.2 %). This is our estimate of Mp. 

The combined uncertainty in the refractive index, u„, 
is obtained as the quadrature sum of the standard devi- 
ation of the mean and the two Type B uncertainties. The 
resulting value is 0.0061. The effect of this refractive 
index uncertainty on the uncertainty in the diameter of 
SRM® 1690, Mbi, is determined to be 0.0020 ^m based 
on the analysis on page 14 of the study by Mulholland 
etal. [1]. 

The combined uncertainty for the mean diameter is 
computed using Eq. (4) with the corrected values of Mgb 
Mb4, and Mb5- The resulting value is 0.0026 |im. Equation 
(5) is used for computing the degrees of freedom with 
a resulting value of 1.47 x 10^. For a 95 % confidence 
level, the corresponding coverage factor is 2.00. 
Therefore the corrected expanded uncertainty is 0.0052 
|im. This value is about 1/3 less than the value of 0.008 
^m on the SRM® 1690 certificate. 



second Type B component, m„ is the uncertainty in the 
determination of the point in the particle image that 
corresponds to the actual edge of the particle. The esti- 
mated value [4] is 0.001 |im. 

The combined uncertainty, obtained from the quad- 
rature sum of Ma> "m> and u^, is equal to 0.00184 |im. 
The effective number of degrees of freedom computed 
using Eq. (5) is equal to 14. In this case the coverage 
factor is 2.20 and the expanded uncertainty is equal to 
0.0040 iim. This value is about 40 % smaller than the 
value currently given on the SRM® 1691 Certificate. 



4. Conclusions 

• The revised expanded uncertainty (approximately 
95 % confidence level) for SRM® 1690 is equal to 
0.005 |im with number mean diameter of 0.895 |im 
compared to 0.008 |im on the SRM® certificates 
dated 2004 and earlier 

• The revised expanded uncertainty for SRM® 1 69 1 is 
equal to 0.004 |im with number mean diameter of 
0.269 iim compared to 0.007 ^m on SRM® certifi- 
cates dated 2004 and earlier 



3. Impact on the Certified Values for 
SRM® 1691 

The uncertainty of the 0.3 |im SRM® is recomputed 
to include the effect of the change in the uncertainty in 
the 1.0 nm SRM®. The current NIST Guidelines for 
expressing uncertainty [2] are used in carrying out the 
analysis. The particle sizes for the 0.3 |im spheres were 
measured by transmission electron microscopy (TEM). 
Both the 1.0 ^m SRM® particles and the 0.3 ^m parti- 
cles were deposited on five TEM grids. For each grid, 
at least 40 of both the 0.3 |im particles and of the 1.0 
|im particles were sized. The 1.0 |im SRM® served as 
the magnification standard for the measurements. For 
each of the five TEM grids a mean size was computed. 
The average of these five mean sizes was found to be 
0.269 |im. The standard deviation of the means was 
found to be 0.00134 |im, which is equal to the Type A 
uncertainty of the measurements, %. 

One Type B uncertainty component is the uncertain- 
ty in the magnification. As shown in the study by 
Lettieri and Hembree [4], the magnification uncertain- 
ty is equal to the combined standard uncertainty in the 
1.0 lim SRM®, 0.0026 ^m, multiplied times the ratio of 
the diameter of the 0.3 ^m SRM® to the 1.0 ^m SRM®. 
The resulting value of Mj^ is equal to 0.00078 |im. The 
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